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rectangular  and  polar  coordinates,  respectively. 

x+iy^re1®,  complex  variable. 

components  of  the  stress  and  strain  tensors, 
respectively . 

components  of  the  surface  traction,  displacement 
and  body  force  vectors,  respectively. 

rectangular  stress  components. 

rectangular  displacement  components. 

crack  dimension. 

crack  inclination. 

applied  load  biaxiality  ratio,  t  (<»)/t  (»)  . 

xx  yy 

uniform  stress  applied  to  the  outer  boundary  surface, 
elastic  shear  modulus. 

Young's  elastic  modulus. 

Poisson's  ratio. 

(3-v)/(l+v)  for  plane  stress,  (3-4v)  for  plane  strain 
surface  density  per  unit  area. 

elastic  strain  energy  density  per  unit  thickness. 

work  of  forces  applied  to  the  body. 

the  potential  energy  of  the  system  and  the  elastic 
potential  energy,  respectively. 

surface  energy. 

bound  region  of  the  x-y  plane, 
closed  boundary  curves  of  R. 

holomorphic  functions  of  the  conrolex  variable  z. 


ON  THE  FRACTURE  STRESS  FOR  THE  INCLINED 
CRACK  UNDER  BIAXIAL  LOAD 


1.  Introduction 

The  significance  of  load  biaxiality  for  fracture  was  investigated 
initially  for  the  horizontal  crack  [1],  and  then  for  the  inclined  crack 
[2].  Attention  in  each  of  these  studies  was  confined  to  analysis  of  the 
local  region  enclosing  the  crack-tip.  Calculation  of  the  critical  (frac¬ 
ture)  stress  under  biaxial  applied  load,  in  which  the  Griffith  crack 
instability  hypothesis  was  utilized,  then  followed  for  the  horizontal 
single  crack  and  for  horizontally  oriented  collinear  cracks  [3,4].  A 
similar  type  of  calculation  was  also  performed  for  the  so-called  'shear 
panel '  [ 5 ] . 

In  this  paper  the  value  of  the  tensile  vertical  load  necessary  to 
fracture  a  biaxially  loaded  plane  infinite  sheet  (in  plane  stress  or  plane 
strain)  with  an  inclined  crack  is  determined  using  the  Griffith  fracture 
criterion.  The  infinite  plane  geometry  is  used  initially  as  a  prototype 
enabling  analytical  demonstration  of  the  qualitative  features  of  the 
problem  as,  for  example,  the  interplay  of  the  dependence  of  the  critical 
tensile  load,  o^,  on  the  load  biaxiality,  the  crack  orientation  and  the 
Poisson  ratio  for  the  material. 

A  similar  calculation  for  a for  a  plane  body  having  finite  dimensions 
and  an  inclined  central  crack  using  the  finite  element  method,  is  to 
follow  for  comparative  purposes. 

Experimental  data  from  several  biaxial  load  fracture  test  programs, 

[6-9],  using  horizontally  oriented  cracks,  encompassing  three  different 

kinds  of  structural  materials,  are  summarized  here.  This  data  provides 

some  measure  of  corroboration  of  the  qualitative  features  of  the  dependence 

of  the  critical  tensile  load  upon  the  load  biaxiality  and  the  material 

Poisson  ratio  as  indicated  by  the  prototype  analysis. 

Manuscript  approved  October  3,  1986. 
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2.  Analytical  Solution 

The  stress  boundary  conditions  of  the  plane  infinite  cracked-body  shown 
in  Fig.  1  are,  relative  to  the  (x'-y')  coordinate  system. 


t '  =*  a  ,  t '  =  ka  ,  t'  =  0 

yy  xx  xy 


as  z  -►  ® 


(2.1) 


whereas  relative  to  the  (x-y)  coordinate  system  they  become 


t  *  ij(l+k)  --Lj(l-k)  cos2a 
XX  2  2 


t  =*  4a(l+k)  +-~j(l-k)  cos2a 
yy  2  2 


as  z  -*■  < 


(2.2) 


t  ■  -^cr(l-k)  sin2ct  . 
xy  2 


Along  the  interior  crack  border 


Cyy(x,0)  *  txy(x’0)  “  0  •  M  <  I a|  •  (2.3) 

As  a  problem  of  plane  linear  elasticity,  the  formal  unique  solution 
of  the  boundary  value  problem  posed  by  Fig.  1  is  expressed  in  terms  of  the 
following  sectionally  holomorphic  functions  of  the  complex  variable 
z  »  x  +  iy  =■  re19  [10], 


where 


$(z) 


2  2.J, 


( 2  — a  ) 


iHz) 


Cl2 

2  2.*  +  C2 


(2  -a4*) 


(2.4) 


Cx  -  |[(l-e2ia)+k(l+e2ia)] -|(l+k)  -2-(l-k)[cos2a+isin2a] 
C_  -  -  ■|'(l-k)e2^a  ■  -  -2-(l-k)  [cos2a  +i  sin2a]  . 


(2.5) 


Needed  in  the  sequel  are  the  integrals  of  i>(z)  and  n(z)  ,  (expressed  to 
within  negligible  arbitrary  complex  constants  that  contribute  only  to  rigid 
body  displacements) , 


*(z) 


r  221, 

<t  (z)  dz  =*  C^(z  -a  )  -  Z^z 


oj(z) 


Q  (z)dz  *  C1(z2-a2)lj  +  Z^z 


(2.6) 


Also  needed  is  the  notational  convention  for  conjugation  of  arguments  of 
functions,  values  of  functions  and,  also,  of  their  derivatives.  Thus  if 
z  =  x  -  iy  and 


0(z) 

*  f1(x,y) 

+  if2 

(x,-y) 
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’p(z) 

=*  fx(x,-y)  +if 

2(x»-y) 

o(z) 

■  f]L(x,y) 

-if2 

(x,y) 

T(z) 

*  Kz)  - 

fr(x. 

-y)  -if 2 

(x,-y) 

rcz) 

■ 

♦ 

P’(z)  - 

do,— v 

S!U) 

Viz) 

dO,  . 

*  di(z> 

y 

V ' (z)  = 

do . . 

r?U) 

The  stress  and  displacement  components  are  determined  from  the  solution 
expressions  (2.4) -(2.6)  by  set  of  relations 

Cxx  +  tyy  *  +  Kz)  ]  -  4Re [<J (z)  ] 

txx'itxy  a  «(2)+!l(7)  +  Cz"I)FW  (2.8) 

2u(u  +iu  )  -  <$( z)  -  w(z)  -  (z-z)<t>(z) 
x  y 
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It  will  be  more  convenient  to  have  expressions  (2.8)  in  the  following  un¬ 
coupled  form: 

t  *  Re[$(z)  ]  +2y  Im[<5>  '  (z)  ]  +  Re[ft(z)  ] 

t  *  3Re[i>(z)  ]  -2y  Im[  i  '  (z)  ]  -  Re[Q(z)  ] 

xx 

t  =*  -2yRe[$ ' (z) ]  + Im[Q(z) ]  -  Im[$(z) ]  (2.9) 

xy 

2yu  =  <Re[<j(z)  ]  -  Re[ui(z)  ]  -2y  ImC'Kz)  ] 

X 

2yu^  =*  jclm[ ( 2)  ]  -  Im[oi(z)  ]  -2y  Re[$(z)  ] 

Expressions  (2.4)-(2.6),  together  with  equations  (2.9),  are  sufficient  to 
determine  uniquely  the  stress  and  displacement  fields  at  every  interior 
point  of  the  biaxially  loaded  infinite  cracked-body  having  the  geometry 
illustrated  by  Fig.  1. 


*.  %  *,  %  *, 
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3.  Power  Series  Representation  of  the  Stress  and  Displacement  Fields 

Needed  for  the  subsequent  calculation  of  the  elastic  strain  energy  of 
the  body  are  the  power  series  representations  for  the  stress  and  displace¬ 
ment  components.  Consequently  the  analytic  functions  $  (z) ,  $’ (z)  ,  Q(z)  , 

~.{z)  ,  4>(z)  and  uj(z)  appearing  in  the  uncoupled  stress-displacement  equations 
(2.9)  are  expanded  about  the  origin  as  power  (Taylor)  series  as  follows. 

For  (a/  j  z  | )  <1, 


<*>  -  (crc2>+ci)T  (if  +l  (f / V  (t 


>■<*>  = -MV!  VtV 

z  z  z 


«*>  -  (f )  <-!(!)  +-M=) + 

a(2)  -  (c^+^jidZ+Kf)4*  ^(f)6n 


(3.1) 


*(z) 


oj(z) 


(CL-C2)z 


(C^+C^z 


■  C1  IFf  161—5/ 

(  ' Z  '  \z  /  \z  / 


After  conversion  into  the  polar  coordinates  (r,3),  the  series  (3.1)  when 
substituted  into  Eqs.  (2.9)  produce  the  following  series  expressions  for 
the  stress  and  displacements  in  the  large,  i.e.,  for  all  points  (r,3)  for 
which  r  >  a. 


5 


i  2 

t  =  2CR  +  j^^  [C^(cos29  +  2sinfisin30)-C^(2sin0cos36)  ] 

4  _  -  6  ) 

+  [C^O^-  cos^®  +  3sin0sin50)  - C^(3sin0cos50) ]  +C(~j  ; 

3  3 ( C-^— C^)  R  —  C Cj_‘+^2^  ^  | C^(cos29  “  2sin6sin39)+C*(2sin20+2sin9cos30)  ] 

4  ^  6)  (3*2) 

+  [CR(-|-  cos49-3sin0sin50)+C^(-|-  sin40  +  3sin9cos56)  ]  +C^~^  f 

(  2 

»  -2C^+j^^  [C^(2sin0sin36  -  cos29)  + CR(2sin6cos36) ] 

+  [C^(3sin0sin50  cos  40)  + CR(3sin0cos50) ]  +C(“^  |  , 


2uux  =  cos9  +  sin0sin20^-  C2  sine  +  sin0cos29^>] 

i-  /a  V/^  /  (1-*) - .  3  _  r^~  /  _( l+K )  ,•  1a  ,  3 


(3.3) 


+  cos30 +•£  sin0sin40)>-  sin39  +-|-  sin9cos40-] 

+C^~^|^~  |  (<-l)  (1+k)  rcos6  +  2(l-k)  r(cos9cos2a  +  sin8sin2ct)  j 

(3.3) 

2uuy  *  j('T’^Ci  *\^2  ~~  cosS  “  3indsia29y>  +  CR<^— ^  sine  -  sin0cos20^>] 

+  ^1  \^8<m~  cos3®  sin8sin4^  +  "  sin38  sin0cos40^>  ] 

+C^~y^j+ ■£■  (1+k)  (k-1)  rsinQ +y  (1-k)  [rcos0sin2ct  -  rsin6cos2a]  . 

In  these  equations 

CR  *  ^  [ (1+k) -(1-k) cos2a] 

C?  *  -  (l-k)sin2a 

(3.4) 

(CrC2)R  *  |  (1+k)  ;  (Cx+C2)R  -  |  (1+k)  -|  (1-k) cos 2a 

3(CrC2)R"  (C14C2)R  •  I  (1+k)  +f  (l-Wcos2a  . 


*»  _  •  ..  •  »  r  *  r.  *  f  *  *,  a  f.  *V  - A_*  ..Vu  '-a  M  -*  -  ^  W"_  ^  jc  ^  if« 


4. 


Griffith  Crack  Instability  Criterion 


The  criterion  for  crack  instability,  or  fracture,  introduced  by 
Griffith  [  9  ,11]  requires  that  at  the  onset  of  fracture,  meaning  here 
rapid  crack  propagation,  the  potential  energy  of  the  system,  P,  must 
acquire  a  stationary  (maximum)  value.  The  necessary  condition  that  this 
be  so  is 


(4.1) 


For  a  body  with  elastic  material  behavior,  the  total  potential  energy  of 
the  system  consists  of  the  crack  surface  energy,  F,  and  the  elastic  poten¬ 
tial  energy,  V.  The  latter  is  the  difference  between  the  elastic  strain 
energy,  U,  and  the  work  (energy)  of  the  applied  boundary  tractions  and 
body  force  in  bringing  the  body  to  its  deformed  equilibrium  configuration. 
Thus 

P  -  V  +  r  ,  (4.2) 


wnere 


\  t  e  dV  - 
2  jk  jk 


dA  + 


b,  u.  dV 
k  k 


j ,k*l,2 ,3. 


(4.3) 


The  T,  are  the  traction  components  applied  over  the  surface  S  bounding 
the  region  R  of  the  body,  which  is  also  subject  to  the  body  force  compon¬ 
ents  b,^.  The  u.^  are  the  displacement  components.  The  linear  elastic  strain 
energy  of  the  deformed  body  is 


U 


(4.4) 


where  t,,  and  e,, 
jk  jk 


are  the  stress  and  infinitesimal  strain  components. 


respectively. 


ft  r*  ft* 


An  application  of  Clapeyron's  theorem  [12],  which  states  that 


j  W*  +  I  bk"kdv  - 2U  ■ 


(4.5) 


allows  the  potential  energy  of  the  system  to  be  expressed  equivalently  as 


p  =  r  -  u  . 


The  condition  for  onset  of  fracture  thus  becomes 


(4.6) 


s  <r-°>  - 0  • 


(4.7) 


Since  Griffith  assumes  that  the  crack  surface  energy  is  equal  to  the  crack, 
surface  area  multiplied  by  a  constant  surface  energy  density,  y»  applica¬ 
tion  of  the  fracture  criterion  (4.7)  reduces  to  a  calculation  of  the  deri¬ 
vative  of  the  elastic  strain  energy  of  the  entire  body. 

At  the  risk  of  appearing  somewhat  pedantic,  the  author  believes  that  it 
might  prove  useful  at  this  juncture  to  offer  several  observations  with  respect 
to  the  fracture  criterion  (4,1)  or  (4.7).  Application  of  the  Griffith  cri¬ 
terion  does  not  in  any  way  imply  the  subsequent  direction  of  crack  propa¬ 
gation,  be  it  at  some  angle  with  respect  to  the  plane  of  the  original 
crack,  or  along  the  direction  of  the  original  crack.  There  is  nothing 
within  the  criterion  that  can  allow  for  this.  The  operation  of  performing 
the  elastic  strain  energy  derivative  dU/aa,  by  the  definition  of  derivative, 
requires  the  use  of  two  incrementally  differing  crack  sizes,  for  the  same 
crack  geometry  and  loading,  in  the  limit  as  the  increment  goes  to  zero. 

This  operation  does  not  suggest  that  when  fracture  ensues  the  subsequent  crack 
extension  must  be  along  the  original  direction  of  the  plane  of  the  crack. 

What  it  does  provide  is  the  value  of  the  rate  of  change  of  the  elastic 
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strain  energy  of  the  entire  body  with  respect  to  change  of  crack  1  ngth 
when  the  crack  has  its  original  size.  Nothing  more  is  stated  o.-  implied. 

Another  point  of  interest  deserves  comment.  Were  the  elastic  energy 
derivative  to  have  the  value  zero,  then  the  fracture  criterion,  as  expressed 
in  the  form  given  by  Eq.  (4.7),  would  lose  its  applicability.  Is  it 
physically  possible  for  dU/da  to  have  the  value  zero?  We  know  from  com¬ 
pliance  type  calculations  that  the  elastic  strain  energy  derivative  can  be 
positive  valued  or  negative  valued  under  'dead  load’  or  'fixed  grips' 
boundary  constraints,  respectively.  If,  furthermore,  the  derivative  itself 
is  a  continuous  function  of  the  crack  length,  then  it  must  also  have  the 
possibility  of  having  a  zero  value  as  well.  It  is  physically  plausible 
to  imagine  that  this  can  occur  if  the  rate  of  the  elastic  strain  energy 
increase  with  change  of  crack  size  associated  with  the  deformation  of  the 
body  near  the  outer  boundary,  is  just  equal  to  the  rate  of  elastic  strain 
energy  decrease  (release)  associated  with  the  relaxation  of  strains 


around  the  crack. 


m 

yai 

<y».l 

v'.'.l 


Other  than,  being  simple  (non-self  intersecting)  nested  closed  curves  that 


enclose  the  crack,  A^  and  are  arbitrary  otherwise.  For  convenience 
they  may  be  chosen  such  that  A^  is  a  circle  of  radius  rg>|a|  centered  at 
the  origin  of  the  coordinates,  while  A2  coincides  with  the  borders  of  the 
crack  [cf.  Fig.  3].  Since  line  integration  along  the  crack  borders  re¬ 
quire  that  dy  be  zero,  while  the  boundary  condition  (2.3)  must  also  hold, 
we  have  immediately  that 

I (A.)  -  0  .  (5. A) 


The  line  integral  about  the  outer  circular  boundary  can  be  written  in  the 
form 


I(A1)  *  2  r0 


([t  u  +t  u  ]  cos9  +  [t  u  + 1  u  ]  sin9)d0  . 
XX  X  xy  y  rQ  yy  y  xy  x  rQ 


-IT 


(5.5) 


The  elastic  strain  energy  of  the  infinite  plane  body  can  therefore  be 
expressed  as  the  limit 

U(-)  -  ^  1(A)  +  I  (A  )  -  I  (A.)  .  (5.6) 

0  0 

Substitution  of  the  series  representations  Eqs.  (3.2)  and  (3.3)  for  the 
stress  and  displacements  into  line  integral  (5.5)  leads  to  the  following 
expression 

I(A1}  “I  l6t\  [^(Sjk.a.Kj+a^eja.^+a2  )  Hn(9  ;k,a,<)  ]  d© 

'ir  (5.7) 

where 
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(5.8) 


F(9;k,a,*:)  *  A^[A^cos  Q+A^Ccos  0cos2a  + sin0cos0sin2a) ] 

2  2 
+  A- [A, sin  0 + A,  (sin0cos0sin2a  -  sin  0cos2a) ] 

2  4  5 

2 

+  A^LA^sinOcosa  + A^Ccos  9sin2a  -  sin0cos9cos2a) ] 

2 

+  A^fA^sinacosa + A2(cos0sin0cos2a + sin  0sin2a) ]  , 

G(9;k,a,<)  »  A^(A2 [A^cos^a  +  sin0cos0sin29 ]+A2[AgSin9cos9+sin0cos0cos29 ]) 

+  A2(A2[AgSin^9  -  sin^9cos29 J-A^tA^sinacosa  -  sin^9sin29]) 

2 

+  A^(A2 [AgSin0cos0  -  sin0cos0cos29 l-A^fA^cos  a-sin0cos0sin20 ]) 

2  2  2 

+  A2(A2[AySin0cos9  +  sin  asi^Oj+A^tAgSin  9  +  sin  9cos29]) 

+  y[A^cos0  + A^(rnfi9rns2a-fsin9sln2a) ] (A2[cos0cos2S-2sin0cos0sin30 ] 

-  A^( 2cos9sin29  +  2sin0cos0cos30]) 

1  2 
+  — (A^sin0  +  A^Ccosasii^a  -  sin0cos2a)  ]  (A2[sin9cos29  +  2sin  9sin39  ] 

+ A2[2sin^9cos30])  (5.9) 

+■  y{A4sin0  +  A^cosflsinla  -  sin0cos2a) ] (A2(2sin0cos0cos30 ] 

-  A2[2sin9cos0sin30  -cos9cos20]) 

1  2 
+  y[A^cos9  +  A2(cos0cos2a  +  sin9sin2o)  l^^sin  9cos30] 

-  A2[2sin^9sin30  -  sin0cos20])  , 


A^  »  (l+k)+(l-k) cos2a  ,  A2  3  (l+k)-(l-k)cos2a 


A-  *  (l-k)sin2a 


,  A,*  (k-1) (1+k) 


*7 


jj'  The  functions  H  (0;k,a,K)  appearing  in  the  integrand  of  (5.7)  are  com- 

E 

u  posed  of  sums  and  products  of  sine  and  cosine  functions  and  have  bound 


It  follows  from  Eqs.  (5.6),  (5.11)  and  (5.12),  that  the  elastic  strain 


energy  of  an  infinite  plane  body  which  has  an  inclined  central  crack  and 
which  is  load  biaxially  along  its  remote  outer  boundaries,  can  be  expressed 
as  the  sum 


U(-) 


liB  KA,) 

r0~  1 


UQ  +  AU 


(5.13) 


where 


U, 


lim  r 
ro~  L 


2  a  tt(H-v)  /,  2 


8E 


{* 


k  (l-hc)+2k(ic-3)+<+l 


)] 


(5.14) 


represents  the  unbound  elastic  strain  energy  of  the  infinite  body  not  having 
a  crack,  and 


2  2 

AU  -  T3e~H,)  a-  [k2+2k(2-«)+l-(l-k2)cos2a] 


(5.15) 


is  the  change  of  the  elastic  strain  energy  associated  with  insertion  of  a 
crack  of  length  2a  inclined  at  an  angle  a  relative  to  direction  of  the 
tensile  loading  a. 

When  the  crack  is  horizontally  oriented, a »  tt/2,  Eq.  (5.15)  reduces  to 
the  previously  obtained  form  [  3  j 


AU 


2  2 
q  rr(l-H>)a 


(-i) 


4E 


[ l+k(2-<) ] 


(5.16) 


For  the  crack  with  a  vertical  orientation,  a  =  0,  and  biaxial  loading,  the 
form  for  AU  should  be  identical  to  Eq.  (5.16)  because  of  the  symmetry  of 
the  two  situations.  Thus  for  a»0,  letting  (+k)a  =  t  so  that  a  -  x/(+k) 

■  (+k’)x,  substitution  in  Eq.  (5.15)  recovers  the  form  given  by  Eq.  (5.16), 


AU 


7)  *(1+v)' 


(a«  0) 


8E 


[2k  +2k(2-<)] 


2  2 
T  7T(l+v)a 

4E 


[1  +  k’ (2-tc) ]  , 


(5.17) 


where  x  is  now  the  tensile  load  in  the  direction  perpendicular  to  the  plane 
of  the  crack  and  k'x  is  the  load  parallel  to  the  crack. 

It  is  apparent  from  Ea.  (5.15)  that  the  elastic  energy  change  AU  can 
have  the  value  zero  when  the  horizontal  load,  the  crack  orientation  and 
the  Poisson  ratio  are  such  that  the  condition 

k2  +  2k(2-<) +l-(l-k2)cos  2  a  =*  0  (5.18) 

is  satisfied.  There  are  several  combinations  of  values  of  the  parameters 

k,  3,  and  v,  that  will  satisfy  this  condition.  Among  the  more  obvious  are 

v  ,  1 

(  a  *  2  ’  k  =  <-2 

AU  -  0  :  S  (5.19) 

\  a  *  0  ,  k  =  0 

For  the  first  case  the  crack  is  horizontal  and  the  horizontal  load 
parameter  has  values  that  depend  on  the  Poisson  ratio  of  the  material, 
as  shown  in  Table  1. 

TABLE  1.  Values  of  Load  Biaxiality  and  Poisson's  Ratio  for  Which  AU  *  0 
When  the  Crack  is  Horizontal  (a=»-r/2)  . 


Plane  Stress 

Plane  Strain 

Poisson  Ratio  v 

Load  3iaxiality  k 

Poisson  Ratio  v 

Load  Biaxiality  k 

0.25 

5.00 

0.25 

: k 

0.30 

13.00 

0.30 

-5.00 

0.33 

k 

0.35 

-2.50 

0.40 

-7.00 

0.40 

-1.67 

0.45 

-4.14 

0.45 

-1.25 

When  the  Poisson's  ratio  has  the  particular  values  v  =  0.33  for  plane 
stress,  and  v  ■  0.25  for  plane  strain,  the  elastic  strain  energy  change  AU 
associated  with  insertion  of  the  crack  cannot  have  the  value  zero,  regard¬ 
less  of  the  value  of  the  load  biaxiality  k. 

For  the  latter  case,  AU  ■  0  when  k  =  0  and  a  ■  0,  that  is,  when  the 
applied  load  is  uniaxial  tension  a  and  the  crack  is  oriented  parallel  to 
the  direction  of  this  tensile  load.  If  we  look  upon  this  as  the  limit 
situation 

lim 

ar+0  (AU)  -  0  ,  (5.20) 

k-O 

then  for  small  values  of  a  and  k,  representing  crack  inclinations  close 
to  the  vertical  direction  and  horizontal  loads  of  small  magnitudes,  the 
values  for  AU  will  be  correspondingly  small. 

For  any  of  the  circumstances  for  which  the  parameters  k,  a,  y  are 
collectively  such  as  to  satisfy  Eq.  (5.18),  insertion  of  the  crack  (or, 
mathematically  speaking,  a  slit)  into  the  biaxially  loaded  plane  body 
will  not  cause  any  change  of  the  total  elastic  strain  energy  of  that  body. 
Hence,  as  in  Eq.  (5.13)  is  independent  of  the  crack  dimension,  dU/da 
*  d(Ug  +  AU)da  *  0  in  Eq.  (4.7),  causing  the  Griffith  fracture  criterion, 
as  expressed  by  this  form,  to  become  inapplicable. 
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6. 


The  Critical  or  Fracture  Load 


The  critical  tensile  load,  o  ,  is  obtained  by  application  of  the  Grif¬ 
fith  crack  instability  criterion  (4.7),  where  T  *  4ya  is  the  crack  surface 
energy,  and  expressions  (5 .13) — (5 -15)  represent  the  total  elastic  strain 
energy  of  the  infinite  body.  Thus  the  condition  for  the  stationary  value 
of  the  potential  energy  of  the  system 


_d_ 

da 


[r-u(«)  ] 


_d_ 

da 


[4Ya-(UQ+AU)]= 


0  , 


(6.1) 


leads  to  the  following  expression  for  the  critical  tensile  load  for  an 
inclined  crack 


a 

c 


f  16EV 
(  7ra(l+v) 


[k^+2k(2-<)+l-(l-k^) cos2ct] 


(6.2) 


In  addition  to  the  crack  size,  the  values  of  the  elastic  modulus  and  the 
crack  surface  energy  density,  the  fracture  load  according  to  this  criterion 
is  also  governed  by  three  additional  parameters,  namely,  the  horizontal 
load  (or  load  biaxiality) ,  the  crack  orientation  and  the  Poisson  ratio  of 
the  material. 


ac  =*  f(a,a,k,v,E,y)  .  (6.3) 

When  the  parameter  values  are  such  that  Eq.  (5.18)  holds,  as 
given  by  Eq.  (6.2)  appears  to  become  unbounded  in  value.  Actually  this 
corresponds  to  the  condition  that  AU  ■  0,  with  corresponding  inappli¬ 
cability  of  the  fracture  criterion. 


Figures  4  to  17  illustrate  graphs  of  the  expressions 


(1+v)  [k2+2k^gi)+l-Cl-k2)cos2a] 


1/2 


and 


(1+v) [k2+2k(4v-l)+l-(l-k2)cos2a] 


1/2 


(plane  stress)  (6.4) 


(plane  strain)  ,  (6.5) 


for  different  values  of  the  load  biaxiality,  the  crack  orientation  and  the 
Poisson  ratio  of  the  material.  Symmetry  considerations  with  respect  to  the 
vertical  tensile  load  direction  reduce  consideration  of  the  crack  orienta¬ 
tions  to  the  range  0°  <  a  <_  90°. 

When  the  crack  is  horizontally  oriented,  a  =»  tt/2,  the  fracture  stress 

variation  with  load  biaxiality  is  monotonic  for  all  Poisson  ratio  values 

for  both  plane  stress  and  plane  strain,  as  shown  by  Figures  4  and  8.  There 

is  a  qualitative  change  of  the  variation  of  a  with  k  as  the  Poisson  ratio 

c 

value  varies  from  v  »  0.25  to  v  *  0.45.  For  plane  stress,  increases  ar 
k  increases  in  the  tensile  sense  for  values  of  v<0.33,  while  it  decreases 
as  k  increases  in  the  tensile  sense  when v> 0.33.  A  similar  pattern  exists 
for  plane  strain,  except  that  the  dividing  value  for  the  Poisson  ratio  is 
now  v  *  0.25.  Referring  to  Table  1,  for  v  ■  0.33  and  plane  stress,  or  v  =» 
0.25  and  plane  starin,  the  elastic  strain  energy  associated  with  insertion 
of  the  crack,  AU,  will  have  a  non-zero  value  for  all  values  of  k.  This  is 
also  shown  in  Figures  4  and  8.  Table  1  also  lists  the  collective  sets  of 
values  of  k  and  v  at  which  AU  ■  0.  For  these  sets  of  values  the  denomin¬ 
ator  of  Eq.  (6.2)  has  the  value  zero,  giving  the  apparency  of  ac-*®.  This 
is  illustrated,  for  example,  in  Fig.  8  for  v  *  0.45  with  the  load  biaxiality 
approaching  the  value  k  *  -1.25. 


As  Che  crack  inclination  angle  is  decreased  coward  Che  vertical  orien- 
CaCion,  Che  curves  of  Figures  5-7,  9-11  and  12-17  show  a  peak  effecC  devel¬ 
oping  in  the  vicinity  of  uniaxial  loading,  k  =  0,  becoming  increasingly 
pronounced  as  Che  crack  angle  gets  smaller  and  approaches  a  direction  paral¬ 
lel  to  the  tensile  load  a.  When  the  applied  load  is  uniaxial  tension  and 
the  crack  angle  is  parallel  to  the  load  direction,  the  limit  condition 
(5.20)  holds  and,  correspondingly,  the  fracture  load  implying  that 

it  is  not  possible  to  cause  fracture  under  these  circumstances .  It  follows, 
therefore,  that  for  crack  orientations  that  are  close  to  being  vertical 
and  with  load  biaxiality  small,  giving  a  loading  condition  that  is  almost 
uniaxial  tension,  the  quantity  AU  will  be  small.  Correspondingly  the  value 
for  the  fracture  stress,  ,  becomes  comparatively  large  [cf.  Eqs.  (5.15) 
and  (6.2)].  This  peak  effect  appears  to  be  in  accord  with  the  expectation 
that  it  becomes  more  difficult  to  effect  fracture  (i.e.,  larger  values  for 
c  )  when  the  crack  is  aligned  almost  parallel  to  the  tensile  vertical  load 
direction  while  the  horizontal  load  is  zero,  or  close  to  zero.  According 
to  Fig.  17,  this  effect  is  most  pronounced  for  a  material  with  a  high  Poisson 
ratio  value  in  a  condition  of  plane  strain. 

Figures  12-17  illustrate  that  at  the  particular  load  biaxiality  of 
equal  tension-tension,  i.e.,  k  =  1.0,  the  orientation  of  the  crack  is  im¬ 
material  to  the  value  of  the  critical  tensile  load.  This  is  the  case  for 
all  values  of  the  Poisson  ratio  for  both  plane  stress  and  plane  strain. 

The  difference  in  values  of  the  fracture  stress  between  plane  stress 
and  plane  strain  conditions  (sheet  thickness)  are  greatest  for  the  hori¬ 
zontal  crack  orientation  as  the  load  biaxiality  intensifies  in  tension 
for  plane  stress  and  in  compression  for  plane  strain. 


For  each  crack  oriencation,  for  both  plane  stress  and  plane  strain, 
there  is  a  corresponding  horizontal  compressive  load,  ranging  between 
-0.5  <_  k  <  0,  at  which  the  fracture  stress  has  the  same  value  for  all 
values  of  the  Poisson  ratio,  as  shown  in  Figs.  4-11.  To  the  left  of  this 
cross-over  point  increases  as  the  value  of  v  increases,  whereas  to  the 
right  of  this  point  decreases  as  the  value  of  v  increases. 


.  vv> 


7. 


Fracture  Stress  Test  Data  For  the  Horizontally  Oriented  Crack 


1 


Ever  since  Griffith's  first  tests  measuring  the  breaking  stress  of 


cracked  glass  tubes  in  1921,  there  have  been  only  a  few  other  fracture 


test  programs  in  which  the  loads  have  been  applied  biaxially  [6-9]. 


Moreover,  in  all  of  these  tests  the  crack  orientations  were  parallel  to 


the  horizontally  applied  load,  i.e.,  a  =  n/2.  The  fracture  load  test  data 


obtained  from  these  tests,  although  previously  reported,  are  collected 


and  reproduced  here  because  they  provide  some  measure  against  which  to 


judge  qualitatively  the  results  of  the  analysis,  as  indicated  by  Figures 


4  and  8  for  the  horizontal  crack  orientation.  The  tests  were  Derformed 


on  glass,  several  aluminum  alloys  and  several  plastics. 


Glass  [9 ] 


The  test  specimens  consisted  of  cylindrical  glass  tubes  having  a 


0.02  in.  wall  thickness,  3nd  inner  diameters  that  ranged  from  0.59  in. 


to  0.74  in.  The  tubes  were  cracked  along  the  axial  direction,  annealed 


and  broken  by  internal  pressure  and  end  loading.  The  values  of  the  elas¬ 


tic  modulus,  Poisson  ratio  and  surface  energy  density  for  the  glass  were 


measured  (directly  and  indirectly)  to  be :  E  =  9.01x10  psi,  v  =  0.25, 


Y  =  3.1x10  psi.  The  load  biaxiality  that  was  attained  in  these  tests 


ranged  only  from  k  =  -0.9  to  k  =  0.5. 


Along  side  of  the  test  results  shown  in  Fig.  18  are  Griffith's  1921 


calculation  of  the  critical  stress 


(1921) 


and  his  subsequently  revised  calculation  (for  plane  stress) 
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^n  both  calculations  :  is  independent  of  the  load  biaxialitv,  and  in  the 

c 

second  calculation  it  is  also  independent  of  the  Poisson  ratio  of  the 
. £ 

material.  Also  shown  is  a  plot  of  Eq .  (6.2)  over  the  test  range  of  the 
load  biaxialitv,  which  for  a  =  -r/2  and  plane  stress  becomes 

-  _  (  8Ev_  1 

'c  \  "  h(3v-l)-rl+ 

In  the  absence  of  any  reason  to  think  otherwise,  it  would  be  easy 
to  interpret  this  test  data  as  indicating  a  null  effect  for  the  load  biax- 
iality  on  the  critical  stress,  in  agreement  with  Griffith's  calculations, 
qualitatively  at  least.  :.v"nile  it  would  be  difficult  to  draw  any  conclu¬ 
sions  from  the  limited  data  shown  in  Fig.  18,  it  appears  nevertheless  that 
Eq .  (7.3)  comes  closer  to  fitting  this  data  than  do  either  of  the  Griffith 
formulas . 

Aluminum  Alleys  [6,7] 

Plane  stress  fracture  toughness  tests  from  under  biaxial  loading 

conducted  on  7075-T6,  2024  T-3  and  6061-T4  aluminum  alloys.  Specially 

designed  sheet  speciemns  were  used  in  a  loading  rig  that  was  designed  to 

apply  loading  independently  in  two  directions.  The  specimens  were  0.063 

in.  in  thickness  and  had  crack  length  to  width  ration,  ( 2a/W) ,  ranging 

between  0.40  to  0.50.  The  Poisson  ratio  of  these  alloys  is  approximately 

0.30.  The  critical  tensile  stress  c  was  calculated  from  the  reported  K 

c  c 

* 

A  critical  discussion  of  these  calculations  mav  be  found  in  Ref.  [3], 


values  by  means  of  the  expression 


where 

£(ir)  - 11  -  +  (¥)2] ' 1-20 

is  the  finite  width  correction  factor  used  in  the  test  program. 

7075-T6 

The  results  of  two  test  series  for  this  alloy  are  shown  by  Figures 
19  and  20,  where  the  load  biaxiality  ranged  from  k  *  0  to  k  *  1.8. 

Higher  k  values  could  not  be  tested  successfully  because  the  horizontal 
loading  tabs  of  the  sheet  specimens  were  pulled  off  prior  to  unstable 
fracture,  thus  releasing  the  biaxial  constraint.  Both  series  show  a 
steady  increase  of  the  fracture  stress  with  increase  in  the  magnitude 
of  the  horizontal  tensile  loading  up  to  k  *  1.5,  dropping  off  somewhat  at 
the  limit  horizontal  constraint  value  for  tab  pull  out  at  k  *  1.8. 

Successful  biaxial  fracture  testing  at  higher  tensile  constraint  loading 
beyond  the  k  *  1.8  limit  would  require  the  use  of  thicker  sheet  specimens. 

The  difference  in  the  sets  of  values  shown  by  Figures  19  and  20  is  attributed 
to  the  anisotropy  induced  by  the  processing  of  the  material. 

2024T-3 

Figures  21  and  22  show  the  results  of  two  test  series  for  this  alloy. 

For  this  more  ductile  alloy  the  load  biaxiality  was  limited  to  the  maximum 
value  of  k  ■  1.5  because  of  loss  occasioned  by  horizontal  tab  pull  out. 

Both  series  show  an  increase  of  the  fracture  stress  as  k  increases  k  =  0.0 


to  0.5.  The  fracture  load  then  decreases  to  a  value  at  k  =  1.5  approxi¬ 
mately  equal  to  the  value  of  a  at  k  =  0,  i.e.,  at  uniaxial  tensile 

c 

loading. 

6061-T4 

The  one  series  of  test  results  for  this  alloy  is  shown  by  Fig.  23. 

The  maximum  horizontal  constraint  that  could  be  sustained  successfully 
by  this  alloy  at  0.063  in.  sheet  thickness  was  only  k  *  1.25.  The  pat¬ 
tern  of  variation  of  o  versus  k  is  similar  to  that  for  2024-T3,  with  a 
peak  value  for  a  at  k  *  0.5,  followed  by  a  decrease  to  a  level  at  k*1.25 
approximately  equal  to  the  value  for  at  k  ■  0. 

An  explanation  for  the  appearance  of  the  peak  values  for  at  k  » 

0.5  load  biaxiality  for  the  2024-T3  and  6061-T4  tests  could  not  be  given. 
Since  horizontal  tab  failures  (pull  out)  prior  to  fracture  did  occur  at 
the  higher  k  values,  it  was  conjectured  that  possibly  the  horizontal  tabs 
of  the  test  specimens  for  the  more  ductile  alloys  began  to  yield  prior 
to  the  onset  of  unstable  fracture,  thus  causing  partial  loss  of  the  hori¬ 
zontal  constraint.  However  this  could  not  be  verified  by  direct  experi¬ 
mental  observation. 

Plastics  [6,7,8] 

Three  test  series  were  conducted  on  polyvinylchloride,  PVC,  and 
polymethylmethacrylate,  PMMA,  (plexiglass).  Two  were  run  at  The  George 
Washington  University  using  the  same  biaxial  load  set-up  and  specimen 
design  as  for  the  aluminum  series,  while  the  other  was  performed  at  the 
Imperial  College  of  Technology,  London.  The  test  specimen  thicknesses 
were  0.25  in.  and  0.16  in.  The  Poisson  ratio  for  these  plastics  range 


approximately  between  0.40-0.45.  The  0.25  in.  thickness  specimens  gave 
the  appearance  of  plane  strain  fracture  conditions. 

PVC 

The  fracture  test  results  for  this  plastic  is  shown  by  Fig.  24.  The 
fracture  stress  decreases  steadily  as  the  load  biaxiality  is  increased  in 
tension  from  k  *  0  to  k  ■  1.5 

?MMA 

The  two  sets  of  test  results  for  PMMA  are  shown  by  Figures  25  and 
26.  Both  show  a  wide  degree  of  scatter  of  the  data.  However  the  linear 
regression  lines  for  both  sets  show  to  be  decreasing  over  the  range 
of  load  biaxiality  from  k  =»  -0.50  to  k  **  2.75,  with  the  regression  lines 
almost  parallel. 


3.  Conclusions 


The  test  data  that  is  available  for  horizontally  oriented  cracks  shows 
rather  convincingly  that  the  critical  tensile  load  is  influenced  by  the 
presence  of  the  horizontal  load  and,  also,  by  the  Poisson  ratio  of  the 
material.  Furthermore,  the  nature  of  the  dependency  is,  qualitatively 
speaking,  in  agreement  with  the  results  of  the  critical  load  calculation 

illustrated  by  Figures  4  and  8. 

It  is  apparent  from  the  two  term  nature  of  the  Griffith  global 

fracture  criterion,  Eqs.  (4.7)  or  (6.1),  that  the  criterion  loses  its 
applicability  for  those  circumstances  in  which  the  elastic  strain  energy  of 
the  entire  body  can  remain  unaltered  as  the  crack  is  extended,  i.e.,  dU/da- 
0.  It  is  not  difficult  to  imagine  a  physcially  plausible  situation  where 
this  can  take  place. 

Moreover,  In  application  of  the  fracture  criterion  for  determination 
of  the  fracture  stress  for  an  inclined  crack  in  a  biaxially  loaded  plane 
infinite  body,  calculation  of  the  elastic  strain  energy  of  the  body  also 
reveals  that  there  are  conditions  for  the  load  biaxiality,  crack 
orientation  and  material  Poisson  ratio  at  which  the  derivative  of  the  total 
strain  energy  with  respect  to  crack  size  can  have  the  value  zero.  This 
would  appear  to  indicate  that  the  Griffith  global  fracture  criteria  is  of 
limited  generality,  even  for  entirely  elastic  material  behavior  throughout 
the  body. 
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Fig.  4  —  Fracture  load  versus  load  biaxiality. 


32 


:y)  0  i> 


BIAXIAL  LOAD  PARAMETER  (k) 

Fig.  5  —  Fracture  load  versus  ioad  biaxiaiity. 
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Fig.  10  —  Fracture  load  versus  load  biaxiality. 
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Fig.  13  —  Fracture  load  versus  load  biaxialtty. 
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Fig.  17  —  Fracture  load  versus  load  biaxiality. 
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Fig.  20  —  Fracture  load  for  Alum.  Alloy  707S-T6— plane  stress. 

Crack  perpendicular  to  the  roll  direction.  (George  Washington  Univ.  tests) 
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Fig.  21  —  Fracture  load  for  Alum.  Alloy  2024-T3— plane  stress. 

Crack  perpendicular  to  the  roll  direction.  (George  Washington  Univ.  tests) 
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Fig.  22  —  Fracture  load  for  Alum.  Alloy  2024-T3— plane  stress. 
Crack  parallel  to  the  roll  direction.  (George  Washington  Univ.  tests) 
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Fig.  24  —  Fracture  load  for  PVC— plane  strain. 
(George  Washington  Univ.  tests) 
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Fig.  26  —  Fracture  load  for  PMMA— plane  strain. 
(George  Washington  Univ.  tests) 
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